* / J\ Si * i / . / 

jj b f-G? 5 '//^Y 


f 

Star Formation and the Distribution of HI and Infrared Emission in M51 

Richard J. Rand and Shrinivas R. Kulkarni 
California Institute of Technology 
and 

W. Rice 

Infrared Processing and Analysis Center, 

Jet Propulsion Laboratory, 

California Institute of Technology 


1 Alfred P. Sloan Fellow and Presidential Young Investigator 


J 

(NASA-CR -106928) -> T AR FORMATION AND TrtF N90-2R262 

DISTRIBUTION OF HI AND INFRARED EMISSION IN 
MSI (California Inst, of Tech.) SR p 

CSCL 03A Unci as 

G3/B9 0T026A2 


Abstract 


HI, infrared, CO, Ha and 5-band observations of M51, the prototypical 
grand-design spiral galaxy, are used to study the consequences of star formation 
for the distribution of HI and dust. Using the HI and CO data sets we carry out 
new tests of the idea of Tilanus and Allen that the HI is largely a dissociation 
product in star-forming regions. We confirm that the HI spiral arms are generally 
coincident with the HII region arms, and offset downstream from the CO arms. 
The radial distributions of total gas, Ha and HI surface density have a simple ex- 
planation in the dissociation picture. The distributions also demonstrate how the 
surface density of HI might be related to the star formation efficiency in molecule- 
rich galaxies. The large width of the HI regions along the arms compared to that 
of the HII regions can be understood in terms of a simple “Stromgren sphere” cal- 
culation. The longer lifetime of the stars producing dissociating radiation vs. those 
producing ionizing radiation will also contribute to the greater width of the HI arms 
if stars are continuously forming on the arms. The lack of detailed coincidence of 
the HI and HII regions along the inner arms has a variety of possible explanations 
within the dissociation scenario. 

We carry out two simple tests to probe the origins of the IRAS emission 
in M51. First, we find that the infrared excess (IRE) of M51 is 24, suggesting 
that a substantial fraction of the infrared emission arises from dust outside of star- 
forming regions. Second, radial cuts through the IRAS bands show that at 12, 


- 1 - 



25 and 60 jim, the arm-interarm contrast of the IRAS emission is substantially 

less than that of the Ha emission, convolved to matching resolutions, pro\iding 

further evidence that there are large amounts of dust not associated with star- 
f 

forming regions. Deconvolved IRAS maps have improved resolution but do not 
change the above finding. 

Subject headings: galaxies: interstellar matter - galaxies: individual (M51) - stars: 
formation - infrared: sources 
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I. Introduction 


Great advances in our undemanding of ho. spiral density waves organise 
atomic and molecular gas in galaxies and how they therefore relate to ongoing star 
formation have been made possible in the pas, few years with the advent of high- 
resolution centimeter- and millimeter-wave interferometers. MM has been a popular 
target for such studies because of its proximity (9.6 Mpc; Sandage and Tammann 
,975), nearly face-on aspect, strong spiral structure, and molecule- rich ISM (about 
90% of the ISM is molecular; Scoville and Young 1983). Observations of the ISM of 
this galaxy have been recently reviewed by Rand and Tilanus (1990). 

In two previous papers (Vogel, Kuliarni, and Scoville 1988, hereafter VKS; 
Rend and Kulkami 1990a, hereafter RK) w, have presented mosaic maps of CO 
, - 0 emission from M51 a, 8” resolution made with the Millimeter Interferometer 
of the Owens Valley Radio Observatory (OVRO). We detected molecular spiral 
arms which are well aligned with the dust lanes but offset upstream from the HII 
region arms (see Figure lo.h), implying that the onset of massive star formation 
occurs several million years after the peak compression of the molecular gas. The 
spiral arms represent about onuthird of the total emission measured with single 
dish instruments. The tangential and radial velocity shifts predicted by density 
wave theory as the gas is being compressed can be clearly seen where the arms 
cross the major and minor axes of the galaxy. Using Ho emission as a tracer of the 
massive star formation rate, we found that the star formation efficiency is higher 
on the arms than between the arms, implying .ha, the density wave is triggering 
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excess star formation from the compressed gas. This result has been confirmed in 
single-dish CO observations (Lord and Yonng 1990). The molecular arms contain 
well-organized structures which we termed Gian. Molecular Associations (GMAs). 
The typical mass of a GMA is 3 x 10' M 0 . GMAs are also found between the 
arms. The on-arm GMAs are roughly gravitationally hound, while the interarm 
GMAs are unbound. The CO 1-0 spectra of the GMAs show that they typically 
consist of one to five velocity components. The GMAs may form through collisional 
agglomerations of pre-existing smaller clouds or by gravitational ins.abiU.ies in the 
molecular gas (see Rand and Kulkami 1990b). The interarm GMAs may owe their 
existence to a weak, secondary compression of the density wave. 

VKS and RK examined how the density wave organizes the molecular gas 
and discussed the consequences for star formation from analysis of Ha image 
this paper, we further examine the “back reaction" of star formation on the ISM of 
M51 by considering, in conjunction with molecular data, the distribution of atomic 
gas and infrared emission from two new data sets: the VLA map of 21-cm emission 
by Rot, et of. (1990, hereafter RBHAC), and the IRAS (Infrared Astronomical 
SateUite) pointed observations of M51 (Rice et ol. 1988). 

Recently, evidence has been presented that in the molecule-rich galaxies M51 
(Tilanus and Allen 1989, hereafter TA) and M83 (Allen, Atherton, and Tilanus 1986) 
HI is predominantly a product of dissociation of H, in star-forming regions. If the 
detected HI were distributed in the same way as all the gas we would expect to see 
the peak of the 21-cm emission along the dust lanes, where the highest compression 
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• _ However the observations show that in these galaxies 
of the gas is occurring. However, 

HI is concentrated downstream from the dust lanes but coincident with 

rep0B arms. TA estimate M in Ml ^ ^ “ 

_ i trk t u WSR'l' map of TA, the new 

dissociated to star-forming regions. Compared 

VIA 21-cm emission map of Ml by RBHAC ban higher resolution (S' - ^ * 
18 .) and sensitivity. In this paper we nse this map, along with on. Ho and B-hand 
images, CO data (me RK) and the sto^dish CO data of land and Young (1990) 
lo carry on. further guantitative tests of the dissociation scenario. In particular, we 
examine whether a) the radial distributions of total gas. Ho and HI, b) the small 
scale structure of the HI emission and c, the detailed spatial relationship between 
Urn HI and HII regions can be explained in this picture. We find that the new VLA 
m ap provides further evidence in favor of the dissociation scenario. 

Radiation from yonng stars no. only ionites HI and dissociates Hr but also 

heats the intersteUar dust. The infrared emtosio. from dost allows yet another probe 

.v tcm t„ the second part of this paper, 
of the back reaction of star formation on the ISM. In the s 

, a- „r to Study the distribution and origins 

we use the IRAS pointed observances of M51 to study 

spirals (Cox and Meager 1987; Unsdale-Persson and Helou 1987; Helou 1986) has 
led to a picture for the dust emission in which the 12 pm emission is dominated by 

very small grains transiently heated by the intersteUar radiation held (ISRF; there 

, . , nw /T o stars') the 25 and 60 /im emission is from 

may also be a contribution from OH/IR stars), me 

warm dust in star-forming regions but may also have a substantial contribution 
from generally distributed dust heated by the ISRF, and the 100 pm enuss.on ,s 
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dotted by relatively large grain, heated by the ISRF. The very smaU grains 
„e expected to be destroyed in star-forming repons (Pnget, Leger, and Bonlanger 
,985), so little 12 pm emiasion i. expected from snch region.. Broadly districted 
dost heated by the ISRF has come to he referred to'as “cirrus" (Low et «l. 1984). 

We attempt to interpret the infrared emission of M51 in term, of these two- 
components: dust in star-forming regions and cirrus. Our goal is to understan 
origins of the IRAS emission in M51. In particular, we carry out two test gn 

distinguish between infrared emission arising from star-fonning regions, and emis- 
sion from a cirrus component. We first calculate the “Infrared Excess" (IRE; Merger 
,978), a dimensionless quantity which provides a measure of the fraction of the in- 
frared Snx which arises from star-fonning regions. We then compare cuts through 
the IRAS maps with cuts through Ho and B-band CCD images convolved to the 
IRAS resolutions to determine whether the in&ared emission is distributed in the 
same way as the optical hue emission from star-forming regions. If the infrared 
emission shows a reduced concentration to the spiral arms compared to the Ho 
emission, then one can conclude that a fraction of the infrared emission arises from 
generally distributed dust. We find that the results of both these tests indicate 
that a substantial fraction of the emission detected by IRAS does no, originate in 

stax-forming regions. 

II. HI and Dissociation 


a). Radial Distributions 
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The radial distribution, of the total gas, Ho and HI surface density are shown 
in Figure 2. To estimate the total (Hr and HI) gas surface density, we have used 
the FCRAO map of CO emission at 45” resolution by Lord and Young (1990). 


To convert the CO emission to molecular column density, we assume the Galactic 
conversion factor of « = 3x 10" mol cm- (K km •")-. The FCRAO map is used 
here and below whenever complete information on the molecular gas surface density 
is required. The total gas profile shows a roughly exponential fall-off with radius, 
while the Ho profile shows a similar fall-off, but with a pronounced enhancement 


a. £ 23 150 ”. The HI profile shows a central hole (R 5 30"), a relatively fiat inner 
disk (30” ZR& WO”), a large enhancement in the outer disk peaking at R = 150”, 
and finally a falloff at R > 175”. This profile can be understood in terms of the 
theory of Federman, Glassgold, and Kwan (1979), who consider the shielding of a 
molecular cloud to dissociating radiation by an outer envelope of HI. In this picture, 
the thickness of the dissociated envelope depends on the ratio of the dissociation 
and formation rates of Hr. This ratio has been shown by Hollenbach, Werner and 
Salpeter (1971) to depend on a power of G/n, where G is essentially the flux of 
dissociating radiation times the cross section in the Lyman bands of Hr , and n is the 
density of the gas. Here we are considering the creation of dissociated envelopes in 
star-forming regions by energetic photons from the young stellar associations. Since 
Ho traces the high-mass stam it should also trace the dissociating radiation. The 
fiat HI profile in the inner disk and the large enhancement in HI around R ~ 150 
can then be understood in terms of this simple balance between molecule formation 
and dissociation in star-forming regions. In the inner disk at AS 100”, the Ho 
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and to. *1 gas profiles show a similar faUoff with radius, so that the balaace between 
dissociation and formation changes little with radius and a constant column density 
of HI is produced. At * a 150”, the “excess” of dissociating radiation in the outer 
arms relative to the gas density tips the scales in favor of dissociation, and thus 
more HI is produced. Beyond R = 170”, the Ho enhancement ends, and the HI 
returns to roughly its inner disk value. The radial distributions therefore support 

the dissociation picture. 

An interesting consequence of such an interpretation for the origin of the HI 
is that, since the Ho emission traces the star formation rate, the HI surface density 
is actually related to the star formation efficiency (the star formation rate per uni. 
gas mass). The possibility of using HI as an indicator of star formation efficiency in 
molecule-rich galaxies such as M51 and M83 should be prnsued further. 

b). The Small-scale Structure 

In Figure Ic.d we show the VLA HI map of RBHAC overlaid on our CO map 
and Ho image. It is clear that compared to the CO distribution there is enhanced HI 
emission in the outer (R 1 100”) arms, although the ammuthally averaged column 
density of HI rises to a peak of only 40% of that of Hr as determined from the 
FCRAO CO map (the units of the Hr column density are atoms cm' 1 ). Comparing 
with the OVRO CO map, we find a typical ratio of Hr to HI column density in the 
inner arms of about 16. Therefore the new VLA map confirms the low HI/H, ratio 
on the arms, and the conclusion of TA that only a small fraction of the molecular 

gas is dissociated. 
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From the VLA HI map (Figures lc,d) we note the foUowing: 

i) Over large sections of the inner (. R -S 100”) arms (see F.gure lc) an offset 
can be seen between the centroids of the HI arms and CO arms, hot the broad HI 
emission in man, plac'es extends npstream to the CO arms and often even further 
upstream into the pre-compression zone. 

U) The inner HI arms show a patchy structure and are broader than the Ha 
arms. The typical diameter (FWHM) of the bright HII regions in the inner arms is 
5.5” or 250 pc, while the typical diameter of the HI arms is about 17” or 800 pc, 
although there is much variation in the latter quantity. 

Hi) While parts of the inner arms show a simple morphology in which HI 
envelopes exist around HU regions, there are many HI clouds which have very little 
Ha emission associated with them, as in M83 (Allen, Atherton, and Tilanus 1986). 
Figure 3 is a close-up view of an arm segment which shows the detailed spatial 
correspondence of the CO, 21-cm, and Ha emission. 

iv) There are a few HI clouds in the interarm regions, where there is also 
very little or no Ha emission. It should be pointed out that the total mass of HI 
in the VLA 8”-resolution map is 2 x 10 9 M e , which is only 40% of the mass found 
in the single-dish observations of Appleton, Foster, and Davies (1986) and Rots 
(1980). In fact, RBHAC report that in their lower resolution VLA observations, 
there is broadly distributed HI between the arms at the level of about 10 2 ° cm" 2 
which has been missed in the 8” map. This interarm emission may not make up 


-9 - 



the entire misting Hex, bn. the feet that the total flex is well below .be single-disb 
value indicates that some emission has been resolved on. in .be 8" map. 

We now explore whether these diverse observational results can be explained 
in the framework of the dissociation model. 


i) The ffl-CO Offset 

First, the (act that the HI arms are, like the Ho arms, offset downstream from 
the CO arms supports the idea that the HI is produced in star-forming regions. If 
the HI arms seen in the VLA map simply consisted of gas compressed by the density 
wave, then we wonld expect to see the HI emission coincident with the CO emission. 
However, the fact that in places the HI emission extends npstream from the CO 
emission implies in the dissociation model that a small fraction of the molecular 
gas is dissociated by downstream star formation before it reaches the zone of peak 

compression. 


ii) The HI Arm Width 

We now address the issue of whether the greater width of the HI arms com- 
pared to the Ha arms in the inner disk can be explained by a simple “Stromgren 
sphere” calculation. Certainly the patchiness of the inner HI arms would be ex- 
pected if each emission feature corresponded to a dissociated region around a young 

stellar cluster or clusters. 

In the Appendix, we calculate the expected sizes of the HII and HI regions 
along the inner arms of M51 using the observed Ho luminosity of a typical inner-arm 
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giant HII region to estimate the ionizing photon rate. The effects of dust opacity 
are included in the calculation by using the HH region extinctions published by van 
der Hulst et al. (1988, hereafter HKCR). With our simple model we predict HII 
and HI region diameters of 220 and 800 pc, in good agreement with the observed 
sizes. Therefore, the characteristic size of the HI arms has an explanation in the 

dissociation picture. 

Two perhaps surprising results were found in the calculation presented in the 
Appendix. First, from the extinction study of HKCR, it is difficult to avoid the 
conclusion that most of the extinction occurs in dust between the HII regions and 
the observer, and that therefore there must be large amounts of dust at heights 
greater than 125 pc (the typical giant HII region radius) above the plane of the 
galaxy. Then, for any reasonable dust density and distribution, the gas and dust 
layer in M51 must be substantially thicker than the Galactic layer, having a half- 
width of about 400 pc. Such a large scale-height could be a consequence of stirring 
up of the ISM by the enhanced star formation and (presumably) supernova rate on 
the arms. Higher magnetic and cosmic ray pressures, as a result of star formation 
or directly due to compression by the strong density wave, could also contribute. 
However, such a large value is unprecedented, and is not expected given that the 
velocity dispersion in the HI gas, ~ 10 km s"\ is not as high as one would expect 
for a gas layer with a 400 pc half-thickness. We therefore view this result with 
suspicion. Our calculation is, however, not too sensitive to our inferred scale-height. 
Using a more modest half-thickness of 150 pc, we find an HI region diameter of 600 
pc, still in reasonable agreement with observations. 
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The second surprising result is that we find from two independent methods of 

estimating the dust surface density that the dust-to-ga s ratio (which is a necessary 

number for calculating the H 2 formation rate) is about equal to the Galactic value. 

f 

This near equality of dust-togas ratios in M51 and the Galaxy may be a cause for 
concern given that the typical metallicity in the inner disk of M51 is three times 
the Solar value (Pagel and Edmunds 1981), so that one might expect a higher dust 
abundance. Using a dust-togas ratio of three times the Galactic value gives an HI 

region diameter of 700 pc. 

Our calculation is more sensitive to the gas filling factor than to the above 
two effects. If the filling factor is less than unity, then our calculated sizes will 
decrease. For a filling factor of 0.2, the calculated HI diameter is 585 pc, and for 
0.1, the diameter is 480 pc. Some of the HI emission features axe of this smaller size, 
but it may be difficult to explain some of the larger emission regions if the filling 

factor is low. 

The H 2 reformation timescale, is about 3 X 10 8 <f> yr, where 4> is the 

filling factor. This timescale can be compared to the time between arm passages, 
x - ftp)" 1 , to see whether a steady state between dissociation and reformation 
can be maintained. With the above density wave parameters and rotation curve, the 
time between arm passages is 5 x 10 7 yr at 2 kpc from the nucleus, and 2 x 10 8 yr at 
4 kpc. Thus it is plausible that the dissociated gas can reform into molecules before 
the next arm passage such that a steady state is maintained, provided <f> ~ 0.1. 
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Regardless of variations in the gas parameters discussed above, the important 
result that the HI regions should be substantially larger than the HII regions is 
robust for reasonable parameters. However, the Stromgren calculation may not be 
able to fully explain the characteristic HI arm width, especially if <f> is as low as 0.1. 

Another effect may contribute to the greater HI arm width. Since the hotter 
stars evolve more rapidly, the ionizing radiation from a stellar cluster will disap- 
pear faster than the dissociating radiation, to which B and A stars can contribute 
significantly. Thus, if young stellar clusters are continuously being formed in the 
arm, then the zone over which dissociating stars are found will be broader than the 
zone over which ionizing stars are found. The dissociation zone will also be broader 
than the typical size of a single HI region corresponding to a giant HII region. Some 
of the HI regions, then, may be due to dissociation by multiple clusters formed at 
different times on the arms. The width of the zone can be estimated as follows. 
From simple blackbody approximations, we find that the stars that dominate the 
dissociating flux are late 0 and early B stars, which have main sequence lifetimes 
of tms « 1 - 2 X 10 7 yr. The component of circular velocity perpendicular to the 
arms is approximately (ft rot - Sl P )Rsmi, where t is the pitch angle of the spiral 
arms, and ft P is the pattern speed of the density wave. Tully (1974a) measured t to 
be 18.5°, and determined a pattern speed, scaled to a distance of 9.6 Mpc, of 38 km 
s -i kpc -1 . The rotation curve is essentially flat at v rot = 210 km s (Tully 1974b, 
Scoville and Young 1983). Therefore, at a distance from the nucleus of R = 3 kpc, 
the width perpendicular to the arms of the zone where dissociating stars should be 
found, assuming that they are all formed at about the same spiral phase, is 0.7 ”r MS 
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(Myr) *7-13". Hence the width of the HI arms is likely due in part to this effect. 
The width expected from this effect is larger near the Inner Lindblad Resonance of 
the wave and smaller near corotation. 

Hi) Poor Detailed Spatial Coincidence of HI and HII Regions 

It remains to be explained »by the detailed spatial coincidence of HII and 
HI regions is poor, particularly in the inner regions where the existence of the 
density wave has been established. One possibility, still within the above dissociation 
scenario, is that some of the HI donds on the arms are simply regions of low total 
gas density compared to the typical on-arm density. In terms of the Federman, 
Glassgold, and Kwan (1979) model, for a given Sox of radiation, the HI column 
density is high in these regions because the molecule formation timescale is relatively 
long. In other words, G/n is relatively high. However, a visual inspection of Figures 
lc and Id shows that HI regions not assodated with Ho emission are sometimes 
assodated with CO emission, and sometimes not, indicating that the HI is not 
necessarily found in regions of low total gas density. 

We prefer the following effects as the more likely reasons for of the lack of 
correlation. First, as mentioned above, the ionizing radiation from a stellar cluster 
disappears faster than the dissociating radiation. This effect would not only lead to a 
lack of correlation but would also explain why some HI regions extend preferentially 
downstream from the HII regions (Figure Id). Another possible contributor could be 
the inevitable supernovae, which will clear out the immediate vicinity of the cluster, 
where the Ha emission is concentrated, much more effectively than the regions where 
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the HI is found, much further away from the duster center. The timescale for this 
dealing process may be comparable to the age of ~ 10 7 yr inferred by Brinks and 
Bajaja (1986) for the ~ 500 pc diameter HI holes they found in M31. Finally, if 
some clusters form without many OB stars, ihey may produce HI regions with little 
Ha emission. 


iv) Interarm HI 

The surface density of the smoothly distributed interarm component is 10 20 
cm -2 . Such a column density is roughly that expected for a medium consisting of 
dissodated envelopes of molecular douds, if the envdopes have a filling factor of 
order 0.1-0.2 (see Federman, Glassgold, and Kwan 1979). The interarm HI thus has 
a simple explanation in the dissodation picture. The few interarm douds seen at 
8” resolution do not have much assodated Ha emission. They may be dissodation 
regions around stellar clusters with few OB stars. 

III. Infrared Properties 

The aim of this analysis is to understand whether or not the infrared emission 
in M51 arises purely from star-forming regions. Figure 4 shows the IRAS maps of 
M51, while Table 1 lists some of the global properties of the infrared emission. TV 
and Tc are the resulting temperatures in a two component fit to the IRAS data by 
Rice et al. (1988). Tc is £ 1 a above the mean value of all large optical galaxies 
modelled in this way by Rice et al. (1988). 
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Our first diagnostic is the dimensionless IRE, which is defined as the total 
IRAS flux (t //„ summed over the four IRAS bands, although some authors use the 
60 and 100 pm bands only) divided by the Lyman continuum flux. Mezger, Smith, 
and ChurchweU (1974) calculate the infrared emission expected from an HII region 
for a given ionizing photon rate, and show how this photon rate can be estimated 
from the thermal flux density. The IRE is calculated as follows (Rice et al. 1990): 


IRE = 9.4 (Hz *) 


Firas( 10 13 Wm 2 ) 
S 2 i(mJy) 


( 1 ) 


where F IRAS is the total IRAS flux and S 2 i is the thermal continuum flux at 21 
cm. The IRE is defined such that a value of unity means that all of the infrared 
emission can be explained as arising from dust heated by Ly-c photons which have 
been absorbed and degraded to Lya by gas. An IRE greater than one implies that 
some of the infrared emission arises from dust heated directly by photons from H- 
ionizing stars and from dust heated by stars which provide little ionization. The 
typical IRE for Galactic HII regions is 6 (Myers et al. 1986). For a galaxy, an IRE 
well in excess of the value for Galactic HII regions indicates a substantial amount 
of infrared emission not associated with star-forming regions. 

In order to calculate the IRE, we need an estimate of the thermal contin- 
uum flux at 21 cm. Thermal-nonthermal flux separations have been carried out 
independently by Tilanus et al. (1988) and Klein et al. (1984), with both studies 
finding a thermal flux at 21 cm of about 75 mJy. Thus the IRE of M51 is about 24, 
which is much higher than the typical value for Galactic HII regions, indicating that 
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a large fraction of the infrared luminosity originates outside star-forming regions, 
presumably from the cirrus associated with quiescent molecular clouds and their 
HI envelopes. The only way to escape this conclusion is to invoke a much higher 
dust-to-gas ratio in M51, such that in the HII regions of M51, dust grains compete 
more effectively for ionizing photons than HI atoms relative to Galactic HII regions. 
However, we show in the Appendix that the dust-to-gas ratio is not significantly 
different from the Galactic value. 

Our second test involves comparing the arm-interarm contrast in the IRAS 
maps with that in our Ha image convolved to the resolutions of the IRAS maps. The 
beam of the IRAS detectors is such that the resolution along the in-scan direction 
(P.A. « 45° for M51) is significantly better than that along the cross-scan direction. 
The in-scan resolution of the IRAS observations is sufficient for the measurement 
of an arm-interarm contrast at 12, 25 and 60 ftm. The in-scan and cross-scan 
resolutions are listed in Table 2. If all of the infrared emission arises in star-forming 
regions, we would expect to see the same arm-interarm contrast in the infrared maps 
as in our Ha image. 

Cuts in the in-scan direction of the original 12, 25 and 60 fi m maps are 
presented in Figure 5 along with cuts through Ha and 5-band images which have 
been convolved to the IRAS resolutions using the elongated IRAS beams. Ongoing 
star formation is traced by the Ha emission while the ISRF is best traced by the 
5-band light. The 100 fim cut is not shown because at the low resolution of the 100 
fim data no spiral structure can be seen in either the map or the cut. 
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Features corresponding to the nuclear region and the outer NE and SW spiral 
arms can be seen in all of the cuts in Figure 5. The ann-interarm contrast in the 
12 urn cut is slightly lower than that at 25 /im. It is dear that the ann-interann 
contrast in the IRAS cuts in all three bands is lower than in the Ha cuts but higher 
than in the 5-band cuts. This result implies that in all three bands there is emission 
both from dust in star-forming regions and cirrus. For the warm dust (25 and 60 /i 
m) this is not too surprising given the above discussion of the IRE. This result is 
also understandable for the 12 /im emission if it is dominated by emission from very 
small grains which are destroyed in star-forming regions. 

Deconvolved images of the IRAS pointed observations (see Neugebauer et 
al. 1984) of M51 at all four IRAS bands were produced at the Infrared Processing 
and Analysis Center (IPAC) using a Richardson- Lucy deconvolution algorithm (Au- 
mann, H. H., private communication). We examined maps at various stages of the 
deconvolution procedure and chose ones for the analysis below based on the level of 
spurious structure introduced in the blank regions around M51 by the deconvolution 
process. We judged that maps produced using more than 5 iterations at 12 and 25 
fim, and 20 iterations at 60 and 100 /zm, were possibly unreliable due to this rather 
subjective criterion. The deconvolved maps using the above numbers of iterations 
are shown in Figure 6, and the achieved resolutions are listed in Table 2. Cuts were 
made through these images as before, and are shown in Figure 7, along with Ha and 
5-band cuts at matched resolutions. Increases in the arm-interarm contrast can be 
clearly seen in the 12, 25 and 60 /im maps, but in the 100 /im map, the companion 
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is only partially resolved and the outer arms remain unresolved. The 100 /im cut 
(not shown) shows no evidence for spiral structure. 

The smoothness of the 100 pm map and cut is surprising when compared to 
Smith’s (1982) map at 170 pm with 49” resolution (his Figure 2) made with the 
Kuiper Airborne Observatory. This map, at only slightly better resolution, shows 
clear spiral structure in the NE and SW. The longer wavelength of his observation 
should guarantee a smaller contribution to the emission from dust in star-forming 
regions, since this wavelength traces colder dust, thus reducing the expected contrast 
relative to that at 100 pm. The likely reason for this discrepancy is that the original 
100 pm map had insufficient resolution to detect any spiral structure at all, even in 
the in-scan direction. Naturally, the deconvolution procedure would not bring out 
any spiral structure if there were none in the original data. 

Figure 7 show the same ordering of cuts by arm-interarm contrast as in Figure 
5, and thus point to the same conclusion as that drawn from the original maps: there 
are at least two components to the infrared emission in the 12, 25 and 60 pm bands. 

The deconvolution we have presented is by no means unique. For compari- 
son, we show in Figure 8 a deconvolution of the 60 pm survey data by W. N. Weir 
(personal communication) using the MEMSYS software (based on the Maximum 
Entropy Method). This deconvolution method, unlike the Richardson-Lucy algo- 
rithm, does not assume that the point-spread function is constant over the field. 
The achieved resolution, although it almost certainly varies over the map, is obvi- 
ously much better than that of Figure 6c, and a higher arm-interarm contrast in 
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the outer arms can be seen. Unfortunately, a comparison of this map with the Ha 
image would be perilous since the resolution varies over the map. It is therefore not 
dear how to perform the necessary convolution of the Ho image before the compar- 
ison. Compared to the pointed observations, the IR.AS survey data unfortunately 
has a worse problem with striping in the in-scan direction, and a stripe can be seen 
running through the center of M51 in Figure 8 which confuses the interpretation. 

We emphasize that our condusions which are based on differences m arm- 
interarm contrast between the JJMS and Ha images do not depend on the reliability 
of the Richardson-Lucy deconvolution since the differences can be seen from the cuts 
through the original IRAS maps (Figure 5). 

We can make a very rough estimate of the fractions of infrared emission 
which are associated with star-forming regions and ISRF-heated drrus. Since the 
emissivity per dust grain of the drrus component should be roughly proportional to 
the ISRF intensity, and since the cirrus arises from dust which is associated with the 
general gas distribution, the profile of cirrus emission should resemble the product of 
the H-band profile and the profile of total gas surface density. The emission profile 
of dust in star-forming regions should resemble the Ha profile. For the construction 
of the total gas profile the resolution of the FCRAO CO map is insuffident. We 
therefore use the new map of CO emission at 15” resolution from the Nobeyama 45- 
m and convolve it to the IRAS resolutions. Using this simple model and the profiles 
shown in Figure 5, we estimate that at 60 /im, 75% of the emission comes from the 
drrus component, while 25% of the 25 /xm emission and 75% of the 12 /xm emission 
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is due to cirrus. For M33, by way of comparison, a more careful decomposition 
by Rice el at. (1990) showed that cirrus accounts for J0% of the 12 and 60 pm 
emission, and 10% of the 25 pm emission. For both these then, the 25 pm 

band contains the highest fraction of emission from star-forming regions. 

The results from these two tests imply that a direct interpretation of the 
IRAS luminosity of M51, either the total CMS luminosity or l F m alone, in terms 
of a star formation rate, would be spurious. The same conclusion has been reached 
in studies of many other nearby spirals, where the IRE's range from 13 to 24 (M33 
(Rice el at. 1990]; M101 [Beichman el at. 1987]; and NGC 4565, NGC 891 and 
NGC 5907 [Wainscoat, de Jong, and Wesselius 1987]). Furthermore, studies of M31 
(Walterbos and Schwering 1987) and the Galaxy (Bloemen, Deal, and Thaddeus 
1990) have shown that most of the FIR emission arises from cirrus in these galaxies. 
Finally, Bothun, Lonsdale, and Rice (1989) found for most spiral galaxies in an 
optically selected sample - using feolfm as a diagnostic of dust heating sources 
along with a simple dust-heating model - that dust in star-forming regions and 
cirrus are both important contributors to the total FIR luminosity. 

The position of M51 in the color-color diagram of Helon (1986) confirms the 
above conclusion. In his Figure 2, M51 lies below the Bne where half of the infrared 
e.i.sc derives from star-forming regions, implying that the cirrus contributes the 

majority of the emission. 

These results may not be surprising given that the galaxy is very gas-rich. 
The total mass of gas has been estimated by Scoville and Young (1983) to be 1.2 X 
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10 10 Af o , of which 9 X 10 9 M o is in molecular form. In comparison, Scoville and 
Young estimate a dynamical mass of about 10 n -V/ 3 , so that the gas accounts for 
fully 13% of the dynamical mass. The high SFR in M51 implies significant dust 
heating in star-forming regions, but should also imply a relatively intense ISRF, 
since, for a normal IMF, there will be copious production of intermediate-mass stars 
which dominate the ISRF. This combination of a rich ISM with an intense ISRF 
could be the reason for the significant cirrus component. 

Throughout this analysis, we have used the term ‘'cirrus” to refer to emission 
from dust outside of regions of high-mass star formation. It should be noted, how- 
ever, that Cox, Krtigel, and Mezger (1986) conclude that non-ionizing young B and 
A stars could heat dust in star-forming regions to 30-40 K and provide some of the 
emission detected by IRASin the Galaxy. It has been suggested for our Galaxy (see 
Scoville and Sanders 1988), although there is no direct evidence for M51, that B and 
A star formation is more uniformly distributed than 0 star formation. If this is the 
case in M51, then such an infrared emission component from regions of low-mass 
star formation may have a spatial distribution more like the cirrus component and 
thus be observationally indistinguishable from cirrus based on the spatial distribu- 
tion of far infrared (FIR) emission alone. It is therefore worth keeping in mind that 
radiation from “generally distributed dust” may have more than one origin. 

An important consequence of the main conclusion of this section is that in- 
frared emission cannot be used to independently test the hypothesis of density wave 
triggering of star formation unless the component of emission from star-forming 
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regions can be effectively isolated, both spatially and spectrally. One should keep 

in mind that the arms which have been resolved or partially resolved by IRAS are 

beyond the radius of corotation of the density wave and are probably tidal arms gen- 

f 

erated by the passage of the companion (Tully 1974b; Elmegreen, Elmegreen, and 
Seiden 1989). Therefore, even the deconvolved IRAS images do not begin to resolve 
the density wave arms and hence cannot be used to test density wave triggering of 
star formation. 

IV. Conclusions 

We have used VLA HI data and IRAS maps of M51 in an attempt to further 
understand the consequences of star formation in this molecule-rich galaxy with a 
strong density wave. From the HI map, we find compelling evidence which favors 
the idea that the HI is almost entirely a product of dissociation of H, in star- 
forming regions. The radial distribution of HI as well as the characteristic size of 
the emission features along the arms have simple explanations in the dissociation 
scenario. The lack of detailed correlation between HI and HII regions has a variety 
of possible causes. The low-level interarm HI emission can be explained as being due 
to dissociated envelopes of interarm molecular clouds. There are indications that 
the gas-to-dust ratio is about equal to the Galactic value. An outstanding issue is 
the thickness of the gas and dust layer in M51. From the comparison of measured 
extinctions with models one infers that the layer is substantially thicker than that of 
the Galaxy. However, the large gas dispersions expected in such a case are not seen 
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and there is no direct evidence to support this result. The issue deserves further 
attention. 

The IRE and the reduced arm-interarm contrast of the IRAS maps compared 
to the Ha image both indicate that a substantial fraction of the infrared emission 
detected by IRAS does not arise in star-forming regions, but rather from dust asso- 
ciated with generally distributed gas. Given the high surface density of gas in M51, 
this result is perhaps not too surprising. The same conclusion has now been drawn 
for several nearby galaxies. An important implication of this result is that it will be 
difficult to interpret the FIR emission from M51 in terms of recent star formation, 
and thus to use FIR emission to check the hypothesis of density wave triggering of 
star formation. Regardless of whether a reliable separation of emission components 
can be done, the resolution of the IRAS maps is insufficient to resolve the inner, 

density-wave spiral arms. 
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Appendix 


We show here through a “Stromgren” analysis that the large observed sizes 
of the discrete HI regions along the inner arms can be understood if they owe their 
existence to photo-dissociation by the OB associations along the arms. 

A typical HII region on the inner arms has an extinction- corrected Ha lumi- 
nosity of 5X10 39 erg s' 1 (HKCR), which corresponds to about 3.5x10 s1 ionizing 
photons per second, or 70 05 stars (Kennicutt, Edgar, and Hodge 1989). We as- 
sume that all photons shortward of 13.6 eV are absorbed by H atoms, while photons 
with energies between 11.2 and 13.6 eV are available to dissociate Hj molecules. We 
ignore heavy elements for the moment; however we do allow for dust opacity. 

At any radius r within an HII region, we have the following balance between 
ionization and recombination. 


-^n H ia H ie- (TD(r)+TH,(r)) = ( A1 ) 

4k r 2 

where Q is the ionizing photon rate, n„, is the atomic H density, a H i is the atomic ab- 
sorption coefficient, averaged over frequencies above the ionization threshold, r D (r) 
and fHi(r)) are the frequency- averaged dust and atomic H opacities at radius r, n e 
is the electron density, and a B is the appropriate recombination coefficient in the 
“on-the-spot” approximation. 

For the HI case, we assume that all photons with 11.2 eV < E < 13.6 eV are 
available for the dissociation of H 2 . Using simple blackbody approximations and a 
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Sal peter Initial Mass Function, we estimate that the photon rate in the above energy 
range for a stellar duster is about 50% of the ionising photon rate. We use this ratio 
for the HU regions of MSI. Only about 10% of absorptions of such photons lead 
a dissodation (Spitzer 1978, p. 100), so a factor of 0.1 must be induded in the 
usual Stromgren calculation to account for this. Thus we have 


- - n(Hi)f ? Ua <, ~ (Tp(r)+TH3(r)) = 5n ff n ( HI ) ^ A2 ' 

8(brr 2 

where n(H 2 ) is the molecular density, u„, is the molecular absorption coefficient 
averaged over the above energy range, f„(r) and nt,(D are the dust and Hr opadties 
at radius r averaged over the above energy range, »„ is the total density of H and 
n(HI) is the density of atomic H. R is the Hr formation constant (Spitzer 1978, p. 
123), which is proportional to the projected area of dust grains per H nudeus, so 
that it must be scaled from the Galactic value of Rg,i = 3.0 x Mr 11 cm 3 s ' by 
tbe estimated dust-togas ratio contrast between M51 and our Galaxy (assuming 
similar grain properties). So equation A2 becomes 


J^n(H 2 )a H3 e-( TD(r)+ ^ (r)) = (A3) 

where Dms, and D„„ are the dust-to-gas ratios in MSI and the Galaxy. 

The dectron density can be estimated from the HI1 region study of HKCR. 
They used a 6-cm VLA map with a resolution of 8’, and an He image smoothed 
the same resolution to estimate HII region projected areas and mean densities. 
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Using their measured diameters of about 10” from their smoothed image, they found 
a typical density of 1.5 cm -3 . Since the actual measured diameter of the typical HII 
region in our image is 5.5”, the projected areas determined by HKCR are too large, 
and the density must therefore be scaled up by the correct projected area, giving 
n e = 5 cm -3 . This density is typical of extragalactic giant HII regions (Kennicutt 

1984). 

The extinction in the HI and HU regions can also be estimated from HKCR. 
From a comparison of the extinction at H/? and the Hct-H/1 color excess with the 
models of Caplan and Dehaxveng (1986) for 11 giant HH regions , they deduce 
that only about 0.5 mag of the typical 1.8 mag of visual extinction could arise 
from within the HH regions, and that the remaining extinction occurs in foreground 
dust. A model with most of the extinction occurring in the HII regions appears 
to be strongly ruled out by their data. Assuming spherical symmetry for the HII 
regions, one therefore requires most of the dust to be at heights above ~125 pc, 
which is the radius of a typical giant HII region. It would therefore seem likely that, 
if the dust density above the HII regions is not too different from that inside the 
HII regions, the dust layer half- thickness is perhaps 3-4 times the typical giant 
HII region radius, or ~ 450 pc. Alternatively, one can assume spherical symmetry 
for the HI regions and again infer that the gas, and presumably the dust, extends 
to about 400 pc above the plane. This issue is discussed further in the main part 
of the text since this scale-height is unusually large. We adopt, then, a total Av 
of 0.5 mag in the HII region, and 1.8 mag in the HI region. We use the Galactic 
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extinction curve from Spitzer (1978, p. 158) to convert these to extinctions at 1000 

k. 

The dust-to-gas ratio can simply be estimated from the typical HII region 
visual extinction of about 1.8 mag and the mean gas surface density at 3 kpc from 
the nucleus of 100 M 0 pc' 2 (Lord and Young 1990). This extinction represents the 
contribution from the HII region and the dust layer above it. There should therefore 
be an additional 1.3 mag of extinction from the dust associated with the gas layer 
on the opposite side of the plane, leading to a total extinction of 3.1 mag. Assuming 
the gas and dust are coextensive, these numbers imply a column density of 2.1 X 10 21 
cm' 2 per magnitude of visual extinction, compared to the Galactic value of 2 x 10 21 
cm -2 (Spitzer 1978, p. 156, 161). The inferred dust-to-gas ratio is thus about equal 
to the Galactic value. A check on this number can be obtained from Smith’s (1982) 
estimate of the dust mass surface density in the inner disk of M51 (Smith’s “Bar 
region) of 0.64 M 0 pc" 2 . The dust-to-gas mass ratio calculated in this way is thus 

0.006, again about equal to the Galactic value. 

For n H , we use the above H 2 surface density and find n H = 4 cm" 3 . For 
n(HI), we estimate a value of 0.5 cm" 3 on the arms from the VLA map. Both 
numbers assume the above half-thickness of the gas layer of 450 pc. 

Equations A2 and A3 are easily solved for the ratio of the Stromgren radius 
in the presence of dust to the dust-free Stromgren radius by the method described 
in Spitzer (1978, p. 111-113). Since we used 'the observed sizes to calculate the dust 
opacity, our solution must reproduce these sizes in order to be self-consistent. In 
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fact, we derive diameters of 220 and 800 pc for the HII and HI regions, in good 
agreement with observations. 

Of the most abundant heavy elements, He, N and O will not compete with 

f 

H 2 for the dissociating photons since their ionization potentials are above 13.6 eV. 
Only Cl has an ionization potential (11.3 eV) in the relevant energy range. We 
show here, however, that we do not expect Cl to compete with HI for this spectral 
range of photons. In equilibrium, the absorption rate by H, molecules will equal 
10 times the formation rate, or 3.5 X 10-“ cm' 3 s" 1 . In contrast, essentially all 
absorptions by Cl in this energy range will lead to ionization, and by balancing 
ionizations with recombinations, using a heavy element abundance of 3 times Solar 
(Pagel and Edmunds 1981), we find a photon absorption rate equal to 2.5 x 10 17 
cm -3 s -i Therefore, our results should not be affected by the presence of Cl. 
Including depletion into dust and molecules would strengthen this conclusion. 
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Table 1. In frared Properties of M51 

Reference 


/*( 12 H 

11 ± 5 Jy 

Rice et al. 1988 

/t/(25 /im) 

17 ± 6 Jy 


7^(60 Mm) 

109 ± 4 Jy 

n 

/„(100 Mm) 

292 ± 7 Jy 

n 

/„( 170 Mm) 

393 Jy 

Smith 1982 

Lfjfi (40 Mm - 120 Mm) 

1.9 x 10 10 L , 

Rice et ai 1988 

Lfir/Lb 

0.7 

n 

Liras (5 M m ~l mm) 

5.0 x 10 10 L 3 

n 

Infrared Excess (IRE) 

24 

This work 

Tw 

209.4 K 

Rice et al. 198S 


Tc 


33.4 K 


Rice et al. 1988 



Band 

Original Resolution 
In-Scan/Cross-Scan FWHM 

Deconvolved Resolution 
In- Scan /Cross- Scan FHWM 

12 ;tm 

45”/'271” 

30”/78” 

25 /im 

43”/271” 

32”/81” 

60 /im 

83”/278” 

33” /56” 

100 £tm 

170”/290” 

59” /69” 
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Figure Captions 


FIG. 1. Shown in the top left panel are contours of CO emission from the OVRO 
mosaic overlaid on a red CCD image. The contour levels are as stated in RK. The 
top right panel shows CO contours overlaid on an Ho CCD ,mage. The bottom 
left panel shows contours of CO overlaid on a grey-scale representation of the VLA 
21- cm map. The bottom right panel shows contours of 21-cm emission on the Ha 
image. The contour levels are 1, 1.5, 2, and 4 times 8.4 x lO 10 cm" 2 in the bottom 

right panel. 

FIG. 2. The radial distributions of total gas surface density. Ha emission, and HI 
surface density. 

FIG. 3. A close-up of a spiral arm segment south of the nucleus showing the small 
scale relationship between HI, Ho and CO emission. The solid contours are HI 
emission. The levels are 1, 2, 3, 4, and 5 times 8.4 x 10 ! “ cm" 2 . The dashed 
contours are CO emission from the OVRO mosaic. The grey-scale represents Ho 
emission. The flow direction is indicated. The area labelled “A” is an example of a 
region where an HI peak coincides with an HII region. The areas labelled “B” are 
examples of HI peaks between HII regions. 

FIG. 4. The original maps of IRAS emission in M51: a) 12 M m, b) 25 M m, c) 60 M m, 
and d) lOOjun. The unit contour interval is the estimated noise in the map. 

FIG. 5. Cuts through the original IRAS maps and Ha and H-band CCD images 
along the IRAS in-scan direction from NE to SW. 



FIG. 6. The deconvolved maps of IRAS emission in Mol: a) 12/im, b) 25/im, c) 
60/xm, and d) lOO^m. The unit contour interval is the estimated noise in the map. 

FIG. 7. Cuts through the deconvolved IRAS maps and Ha and 5-band CCD images 

f 

along the IRAS in-scan direction from NE to S\V. 

FIG 8. Deconvolution of the IRAS 60/rm survey data on Mol by W. N. Weir using 
the MEMSYS package. The contour levels are 5, 10, 15, 20, 25, 30, 35, 40, 50, 60, 
70, 80, and 90 percent of the peak. 



Postal Addresses 


Shrinivas R. Kulkarni and Richard J. Rand 
Division of Physics, Math and Astronomy 
California Institute of Technology 
105-24 

Pasadena, CA 91125 
VV. Rice 

Infrared Processing and Analysis Center 
California Institute of Technology 
100-22 


Pasadena, CA 91125 


















NGC519 M51 1 20.HGEOM.1 





quo_i-z<i — oz 


RIGHT ASCENSION 

Peak flux = 8.3567E+00 MJY/SR 
Levs = 1.6559E-01 * ( 3.000, 6.000, 10.00, 

20.00, 40.00) 





NGC519 M51 600.HGEOM.1 


D 

E 

C 

L 

I 

N 

A 

T 

I 

0 

N 



13 28 45 30 15 00 27 45 30 

RIGHT ASCENSION 

Peak flux = 7.571 2E+01 MJY/SR 
Levs - 2.5147E— 01 * ( 3.000, 6.000, 10.00, 
20.00, 40.00, 80.00, 160.0) 





tensity 







zo — i>z-romo 



20 


1 


8 t 


13 28 45 


30 15 00 27 45 30 

RIGHT ASCENSION 
Peak flux = 1.6243E+01 MJY/SR 
Levs * 1.6559E— 01 * ( 3.000, 6.000, 10.00, 
20.00, 40.00) 


1 



zo i>z-rono 


47 36 


NGC519 M51255.HGE0M.1 
II 7P T® I I 


< 



18 

13 28 45 


_L 

30 


15 00 27 45 30 

RIGHT ASCENSION 

Peak flux = 2.8520E+01 MJY/SR 
Levs = 1 .7925E—01 * ( 3.000, 6.000, 10.00, 

20.00, 40.00, 80.00) 


15 


00 


o 







NGC519 M51 1 0020.HGEOM.1 



RIGHT ASCENSION 

Peak flux = 3.0874E+02 MJY/SR 
Levs = 3.6593E— 01 * ( 3.000, 6.000, 10.00 
20.00, 40.00, 80.00, 160.0, 320.0, 640.0) 














